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whose final states are non-CP eigenstates. The penguin contributions to the mixing- 
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induced interference quantities are expected to be about 4% or smaller. We also point 
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out two observables of O(l), which are pure functions of decay constants and form- 
factors, for a clean test of the factorization hypothesis in neutral- P decays into two 
■ heavy charmed mesons. 
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1 Introduction 



In the first-round experiments of P-meson factories, the measurement of Bd — > D + D~, 
D* + D~, D + D*~ and D* + D*~ decay modes will be available. These channels are expected 
to have fairly large branching ratios (of order 10~ 4 to 10 -3 ) and involve the CP-violating 
phase 

where (for % = u,c,t and j = d, s, b) are elements of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix. A determination of (3 from CP asymmetries of Bd — > Z)(*)+Z)(*)~ will be 
useful, not only to cross-check the extraction of (3 from Bd — ► J/tpKs, but also to shed some 
light on the relevant penguin effect and final-state interactions ffl. In addition, it should be 
important to test whether the factorization approximation, which works well for a variety of 
B transitions into one light and one heavy mesons such as and Dp [[J, remains valid for 
the B decays into two heavy charmed mesons. 

In the B®-B® system, the similar decay channels are B s — > DfD~, D* + D~, D+D*~ 
and D* + D*~. The branching ratios of these transitions are expected to be of order 10 -2 , 
therefore they are rapidly accessible at a high-energy B factory. The CP- violating phase 

which is negligibly small in the standard model 0, will be determinable from CP asymmetries 
of B s — > £)(*)+£)(*)"". The measured value of if remarkably larger than 1°, would 
clearly imply the existence of new physics in B®-B® mixing. Whether the factorization 
approximation works well for such decay modes is also an open question. 

As neutral-P decays into CP eigenstates {Bd — + D + D~ , etc.) have been studied in some 
detail |J, the present note concentrates only on Bd — > D*^D T and B s D*^DJ decays, 
whose final states are non-CP eigenstates. These decay modes, in comparison with Bd — > 
P ± 7r T ||, have an obvious advantage for the study of possible CP violation: the amplitudes 
of B° and B° into a common DD* state are comparable in magnitude and therefore large 
interference between them becomes posssible. We show that the weak phases (3 and (3' 
can respectively be extracted from the time-dependent measurements of Bd — > P* ± P =F and 
B s — ► D* decays. A large CP asymmetry is expected to appear in the Bd channels 

2 Note that this result is independent of any specific parametrization of the CKM matrix. Taking into 
account V* b V us + V* b V cs + V t * b V ts = and \V cb V cs \ » \V tb V ts \ > IK&Kal @> one arrives straightforwardly at 

| < axct<m\(V ub V us )/(V cb V cs )\ ~ 1°. 
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within the standard model. The penguin contributions to the mixing-induced interference 
quantities in these decay modes are estimated to be about 4% or smaller. Thus it is a good 
approximation to neglect the relevant penguin effects. We also show that the time-integrated 
measurements allow us to determine two parameters of 0(1), which are pure functions of 
decay constants and form-factors in the factorization approximation: 

t Sd_ A^ D *(mp 

U fn* ' F?* D {ml*) ' 

JDt hi {™d*) 

Comparing the experimental and theoretical results of ( d or ( s will provide a clean test of 
the factorization hypothesis for neutral- 5 decays into two heavy charmed mesons. Finally 
we give some brief comments on possible effects of new physics on B°-B° mixing and CP 
asymmetries. 



2 (3, f3' and CP asymmetries 

The transitions B d — > D* ± D T and B° s — > D^DJ can occur through both tree-level and 
loop-induced (penguin) quark diagrams, The penguin contribution to the overall amplitude 
of each decay mode is tentatively neglected, as its magnitude is small enough in comparison 
with that of the tree-level amplitude (see the next section for a detailed discussion about the 
penguin effect). With the help of this good approximation, one may define two interference 
quantities between decay amplitudes and B d -B d mixing: 

q, A(B° -> D*+D-) V*V td V cb V; d , i{5d+w) 
D * +D ~ ~ P d ' A(B° - D*+D-) ■ ■ V tb V td ' V* b V cd U ■ U 

D ° - q d A(B^D*-D+) V* b V td V cb V: d QdC Qd6 ' W 

where q d /p d = (V t * b V td )/{ V tb V td ) describes the B d -B d mixing phase in the box-diagram ap- 
proximation, ( d and 5 d denote the ratio of the real hadronic matrix elements and the strong 
phase difference between B° d — > D* + D~ and B^ — > D* + D~. For the similar transitions B° s 
vs B° -> D*+D- and D*~Df, we have 

\ = s. A(B° S -» d: + d-) v; b v ts v cb v; s iSs i(5s+2/30 

Dt+D ° - Ps ' A(B° - Dt+D-) - V tb V t * s ' V* b V cs ^ " ^ 

= ^ AjBl -+ D*~D+) V tb V* V; b V cs {Ss _ m 
° rDt ~ <l s ' A(BO - D*~D+) ~ V&V U ' KbK. ^ " ^ ' 1 j 
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in which the notations are self-explanatory. These quantities consist of both weak and 
strong phases and will affect the decay rates in a significant way. Their imaginary parts are 
of particular interest for the study of CP violation: 



ImA £>*+ £,- 
ImA D *- D + 
ImAn*+ 



d: + d: 



ImA 



( d sin(d d + 2(3) , 
( d sm(S d - 2(3) , 
C sm(6 s + 2(3') , 
( s sm(5 s - 2(3') . 



(6) 



It should be noted, however, that ImA/ and ImAj (for / = D* + D~ or D* + D~) themselves 
are not CP- violating observables! Only their difference Im(A/ — A/), which will vanish for 
(3 = (3' = or 7r, measures the CP asymmetry. 

The generic formulas for decay rates of neutral-i? into non-CP eigenstates have been 
given in the literature (see, e.g., Refs. @, |J). For the decay modes under consideration, 
their time-dependent rates read as 



n 



<-> 

B° d (t) 



D* + D- 



oc e 



-r d t 



— ^ + cos(x,r,t) 



( d sm(5 d + 2(3) sm(x d F d t) 



K 



<-> 

B° d (t) 



D*~D~ 



(x e 



-r d t 



1 + CJ (+) 1-CJ 



cos(x d T d t) 



{-) 



+ Cd sin(£ d - 2(3) sai(x d T d t) 



(7) 



and 



K 



(-) 

B° s (t) 



d* + d: 



oc e 



-r,t 



i + C 2 <t> i-C 2 



cos(x s T s t) 



K 



<-) 

5? w 



oc e 



(+) 

- Cs sm(5 s + 2(3') sin(x s T s t) 

Ji-C <+> i-C 2 



cos(x s r s t) 



+ Cd sin(5 s - 2/?') sin(x s r s t) 



(8) 



where x d = 0.723 ±0.032 and x s > 14.0 (at the 95% confidence level) are B° d -B° d and B° s -B° s 
mixing parameters 0, and T d and T s denote decay widths. Note that the four decay rates 
3 Here we neglect the tiny effect from the decay width difference between two B s mass eigenstates. The 
latest theoretical calculation M predicts AT a /T a = 0.054l°;o32, which is negligibly small for our purpose. 



in Eq. (7) or Eq. (8) share a common factor, since the penguin effects in these decays have 
been neglected. 

The time-dependent measurements of — > D* ± D T and B s — ► D^DJ, as shown above, 
allow us to extract the weak phases (3 and (3' up to a four-fold ambiguity: 



sin 2 (2(3) 
sin 2 (2(3') 



1-S+S-) ± J(l- S 2 + )(l- S 2 ) 



(1-S' + S'_) ± J{l-S%)(l-S'l 



(9) 



where S± = sin(5d±2/3) and S'± = sm(S s ±2(3'). Indeed only a two-fold ambiguity in sin(2/3) 
and sin (2/3') remains, because sin(2/3) > and sin(2/3') ~ are expected to be true within 
the standard model Q. The strong phases ^ and 5 S can also be determined from Eqs. (7) 
and (8) up to a four-fold ambiguity: 



sin 2 5 d 



sin 2 5. 



S+S- 
S'S'_ 



± 
± 



1-S 2 + )(1-S 2 ) 



1-S'^)(1-S' 2 ] 



(10) 



If the final-state interactions were insignificant in these decay modes, Sd and 5 S might not 
deviate too much from zero. In this case, S + ~ — £_ and S' + ~ — S'_ would be a good 
approximation. 

The measurement of S± can be carried out at both asymmetric e + e~ 5-meson factories 
and high- luminosity hadron machines. In contrast, a measurement of S± is only available 
at hadronic B factories. With all such time-dependent measurements, the ratios of hadronic 
matix elements Q and ( s are determinable. Note that these two real parameters can acturally 
be determined from the time-integrated measurments at a symmetric e + e~ collider or at 
hadron facilities of B physics. For example, one may integrate the decay rates in Eqs. (7) 
and (8) over t e [0, oo) and obtain ratios of two decay rates as follows: 



K(Bl 



phys 



D* + D~ 



(2 



K(B° 

n(Bl 



d,phys 
phys 



D*+D- 
D*+D: 



K(Bl 



phys 



D*+D: 



« + (2 + x 2 ) 

(2 + xp C 2 + x 2 
x 2 C 2 + (2 + x 2 ) 



11) 



where the subscript "phys" implies that the mixing effect has been taken into account for B® 
and B® mesons (q = d or s). Obviously the first relation of Eq. (11) allows the determination 
of (a- Once x s is known, ( s can be extracted from the second relation of Eq. (11). If a reliable 
calculation of ( s were possible, then the magnitude of x s would reversely be obtained. 
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For the Bd and B s decay modes discussed above, their time-dependent CP asymmetries 
can be defined as 

TZ[B° d (t) -> D* + D~] - n[B° d {t) -> d*-d + ] 



A d (t) 



Tl[B d {t) -> £)*+£)-] + 7e[BS(t) -> D—D- 



Taking Eqs. (7) and (8) into account, we explicitly obtain 

f.s — 2Cd sin(2/3) cos5rf sin(xrfr^) 

dW = (1 + Q) + (1 - CJ) cos(x d r d t) - 20cos(2/5) sin5 d sin(x rf r d t) ' ( } 

and 

, /.x -2Q sin(2g) cos £ s singly) 

sU (1 + C 2 ) + (1 - C 2 ) cos(x a r s t) - 2C cos(2/3')sin<5 s sin(x s r s t) ' 1 } 
Clearly these two formulas would be simplified to the familiar form, if <^( s ) = 1 and 8m s \ = 

held (i.e., if the final states were CP eigenstates). A necessary condition for large CP 

violation is that Q or ( s does not deviate much from one. Subsequently one can see that 

this condition is indeed satisfied - both (d and ( s are of 0(1), estimated in the factorization 

approximation. 



3 Factorization and penguin effect 

Now we calculate the decay amplitudes of B d — > D^D^ and B® — > D^DJ by use of the 
effective weak Hamiltonian and the factorization approximation. The contributions induced 
by the annihilation-type quark diagrams at the tree and penguin levels are expected to have 
significant form-factor suppression and will be neglected. The relevant AB = 1 Hamiltonian 
can then be written as H 



2 10 

Owy EC* Q8 - CWS)£(ciO,) 



+ h.c. 



where q — d or s, q (for % — 1, • 
0{ = (cb) v _ A (qc)v_ A , 
3 = (qb) V -A^2(q'g')v-A , 

5 = (qb) V -A^2(q'q')v+A , 

°7 = 7;(qb) v - A J2[ e q>(<i'q')v+A 



i=l i=3 

10) are the Wilson coefficients, and 



09 = ^(^y-AEhWl 



v-a 



O2 



(g6)v-A(cc)v-A 



loP(i)V-A 



V-A 



(15) 



(g a 6 /3 )v_ A £(g / ( gg^)v + ^ , 

-(g a 6 /3 ) y _ A £ [e/g'^V+A 
<?' 

-(g Q 6 /3 ) y _A£ ^(g^g^y-A] • (16) 
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Here 03,- • -,0q denote the QCD-induced penguin operators, O7, ■ • O10 denote the elec- 
troweak penguin operators, and a and (3 are the SU(2>) color indices. To calculate the 
physical amplitude of an exclusive B decay mode, the Wilson coefficients and the relevant 
hadronic matrix element of four-quark operators need be evaluated in the same renormal- 
ization scheme and at the same energy scale. After this procedure || we are left with the 
effective Wilson coefficients cf , and they will enter the factorized decay amplitude in the 



following combinations [10] : 



S ff C 2^ 1 S S C f 

a 2i = C 2i + -phg J °>2i-l = C 2i-1 + T^ff > ( 17 ) 

c c 

(for i = 1, • • - , 5), where iVf is the effective number of colors. The values of a| ff for differrent 
N^ s can be found in Refs. [10], 111] . 



It is then straightforward to calculate the decay amplitudes of B® — > D* ± D T , B® — > 
D^DJ and their CP-conjugated processes in the factorization scheme. For exmaple, 

A(B° d -> D*-D+) = {(V; b V cd )af - (V t * b V td ) [(of + af ) - (of + of) } M® d + C D *~ , 

- zr+zr) = {(v; b v cd )af - (vsvy (of + a £) } , (is) 

and 

A(5 S ° - £>,-£>+) = {(V; b V cs )af - (V*V ts ) [(of + <) - (of + of) z s ] } M^-- , 

A(fi s ° - D* s + D;) = {(V; b V cs )af - (V t lV ts ) (of + a^)}M5f Dr , (19) 



in which 



Kll = % (m^2)v~A\0)(X*\(bq 3 ) v _ A \B q ) 



V2 

= V2G F m x » f x A* x \m 2 x ) (e x . ■ k Bq ) 
<3 3 = % (X*\(m)v-A\0)(X\(bq s )v-A\B q ) 



V2 

V2G F m x . fx* F* x (m 2 x *) (e x . ■ k Bq ) , (20) 



are the factorized hadronic matrix elements expressed by relevant decay constants and form- 
factors, and 

2ml , , 

z l = 7 TT V 21 

(mi + m 2 ){m b + m 3 ) 

arises from transforming the (V — A)(V + A) currents into the (V — A) (V — A) ones for the 
penguin amplitudes. In Eqs. (20) and (21), ks q denotes the momentum of B q meson (q = d 
or s), and m, is the mass of qi quark (i = 1, 2, 3). 



With the help of Eqs. (18) - (21) we are able to estimate the penguin contributions to 
the interference quantities \d*+d- an d A D .+ D - given in Eqs. (4) and (5), respectively. As 
^> |a^ ff | (for i — 3, • • 10) holds [HI ITTfl, some analytical approximations can reliably 



,eff| 



be made. To the next-to-leading order we obtain 

V tb V td 



IA 



D*+D~ 



1 - 



IA 



d*+d: 



V cb V cd 

v tb v ts 



af + af)z d 

„eff 



COS f3 



v cb v cs 



+ a' 



eflh 



,cff 



cos (5' 



M, 



D+Dt 



Msr D ° 



(22) 



Numerically one finds \(V tb Vt q ) / (V cb V cq )\ ~ 1 for both q = d and q = s cases [||. In addition, 

z s ~ 0.7. Therefore the penguin 



\(af + af)/af 



eft\ 



,eff| 



af/af\ w 5% 0, |TT| and z d 



contributions to \\d*+d-\ an d I^d*+d-I are expected to be about 4% or smaller for arbitrary 
values of (3 and Neglecting the penguin effects in such decay modes turns out to be a 
good approximation, leading to the following instructive result: 



IA 



D*+D- 



|A 



d*+d: 



& 

Cs 



M? d : +D - 

M c D s ! Dr 

mSI +d; 



D 



A^ D *(m 2 



D J 



fD s 



Ft 



B d D, 



III 



D* 



A 



BsD* 




(23) 



As we have emphasized in Eq. (3), ( d and ( s are pure functions of decay constants and form- 
factors. Under SU(3) symmetry ( s = ( d holds. The measurement of these two parameters 
will provide a clean test of the factorization hypothesis. 

At present it remains difficult to make reliable evaluation of the form-factors appearing in 
Eq. (3) or (23). From current data one can determine the ratio /d*//d s = 1.16±0.19 0, but 
little experimental knowledge about fu* / fz> is now available. For self-consistency we just use 



the model-dependent results of decay constants and formfactors in Ref. [12| to illustrate the 
ballpark values of ( d and ( s . We obtain ( d « 1.04 and ( s « 1.03. While these two numbers 
might not be trustworthy, they do reflect the qualitative expectation Q ~ ( s ~ 1. In some 
sense the final states D* ± D =F and D* Df may be treated as "quasi"-CP eigenstates, if the 
strong phases 5 d and 5 S are not significantly large. This feature makes these decay modes 
as interesting as B d -»■ D+D~, D* + D*~ and B s -> D+Dj, D* + D*- for the study of CP 
violation. 

The result ( d ~ 1 implies that the CP asymmetry A d (t) in Eq. (13) can approximately 
be as large as sin(2/3) cos5 d in magnitude. Therefore we are able to get a significant CP- 
violating signal in B d — > D^D^ (e.g., for | cos <5<i| > 0.5) within the standard model. In 
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comparison, the CP asymmetry in B s — > D* Df transitions are vanishingly small in the 
standard model. 



4 Further discussions 



We have discussed CP asymmetries and penguin effects in B d — > D*^D T and B s — > D^DJ 
transitions. A possibility to test the factorization hypothesis for such neutral- P decays has 
also been pointed out. Before ending this work we make two useful comments. 

(1) Beyond the standard model, new physics may introduce an additional CP- violating 
phase into B d -B d or B Q S -B° S mixing. In this case the weak phases to be measured from 
B d — > D*^D T and B s — > D* s Df decay modes will not be (3 and j3'. For illustration we 
only consider the kinds of new physics that do not violate unitarity of the CKM matrix and 



have insignificant effects on the penguin channels of the above-mentioned decays |T3[ . Then 
B d -B d and B®-B® mixing phases can be written as 

Qd = V tb V td i(fap <h = V tb V ts M m (24) 

P d V tb V t * d ' Ps V tb V t * s ' 1 } 

in which NP and 0' NP denote the CP- violating phases induced by new physics. The overall 
weak phases of \d*+d- and A D »+ D - turn out to be 



2(3 



2{3 + (j) 



NP 



2/3' 



2p' + 



(25) 



respectively. Therefore only the combinations {2(3 + 0np) and {2(3' + 4>'^ P ), instead of (3 and 
(3' themselves, can be extracted from the proposed time-dependent measurements. 

In principle it is possible to isolate 0np and from the decay modes Bd — > K°K° and 
B s — > K°K° respectively, which occur only through the penguin- induced quark diagrams. 
Assuming the top-quark dominance, one obtains the relevant interference quantities between 
decay amplitudes and P°-P° mixing as follows: 



Im 



Im 



p d A(B° d ^K°K°) 

■q. A(B° S -> K°K°) 
Ps A(B° S ^K°K°) 



= Im 


f Qd 


v tb v td \ 




{Pd 


V t tV td ) 


= Im 


( <h 


v tb v t * s \ 




\p s 


v; b v ts ) 



sm </> NP 



sin ' 



'np ■ 



(26) 



Unfortunately these decays have quite small branching ratios: B(B d — > K°K°) 
and B(B° S — > K°K°) ~ 10~ 5 |pL 0|| - Hence they could only be meeasured in the second-round 
experiments of P-meson factories. 



Kr b n 
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(2) It is also worth pointing out that the joint decay of a coherent B®B® pair into D* ± Z) =F 
states on the T(4S r ) resonance, i.e., 



1 d,phys d,phys 



T(45) 



(D* + D~ ) Bd (D*~ D + ) Bd 



is interesting, although the decay rate is considerably suppressed (of order 10 8 ). Following 



the general discussions made in Ref. |L4j], we explicitly obtain this joint decay rate (divided 
by \A{B[ 



*+ r>-\\*\ 



D* + D 



as 



n 



T(45) 



OC 



2 + x3 



1 

2 



X ., 



X „ 



\D*+D- 



X 



D—D+ 



X ,.; 



+ 72 - 2 cos(<5 d ) 



1 + ^ 

l + x2 



1 - 



\D*+D~ 



X 



D*~D+ 



sin 2 (2/3) . 



(27) 



We observe that the CP- violating contribution to the decay rate is significant and could 
finally be separated out. Such a measurement might be available in the second-round exper- 
iments of e + e~ I?-meson factories with about 10 9 ~ 10 B®B® events. 

Our conclusion is simply that the decay modes Bd — > D* ± £) =F and B s — > D* Df are very 
useful for the study of CP violation and for the test of factorization. Some special attention 
is worth being paid to them in the upcoming experiments of B physics. 
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